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TR/iN3IENT2 BY GUESTRUGTURINQ '.7ITII DM/iP 

Thomas G. Butlor 
BUTLSR /J^ALYSES 

CUMIVJJY 


Automated Qubatructuring iii level l6 of JIASTRi'JI v;as employed aa a preface 
to the acluticn of o chtroct tranaient analysis. The DMj’iP ALTER atatementc vs:it- 
tcr. to adapt the auboLi-ucturirg for traaiaior.t purpcsec are explained. Data re- 
covery was accomj^liahod v/ith tranafer functiona. Proof o.f the success of the 
method is presoubed vlth on application to a pIgbIc struc-iu’6. 


IMTEODUCTICN 


Bubstructure analysis capability in NABTH.d'I has been automutod for rigid 
formats 1, 2, and 3 only* Patl;er than v.'ait for the oxteusion of automation to 
rigid formc.t 9i it v;as found advantageous to use as much of automated substruc- 
turing lo is now available. Considerable effort v;?s needed tc couple the sub- 
structui'e data to the transient anolyslc and then recover tVie transient responses 
of individual, substructures. 


The tv/o principle features of substructuririg that made this effort woi'th 
v.'Tiile '.;cre the ability to condense small pieces of the mabricus at a time and 
the ability to combine component cuboti'U''tures in different ways. The sbruc- 
tui'o under investigaticn was liable tc damage, and the atbiactic’’ of substruc- 
tui'ing was its ability to substitute damaged oompone;its for the hapless prede- 
ceasorb. Many additional advantages of substructuring became evident during 
the progress of the analysis, but the most notable was the necessity of the 
analyst to organize thoro ugl: J.y , 

The method presented in this paper could be charactorized as a five phased 
£U*,alybis as opposed to the usual three inautomated eubstructuring. Dove Ilerting 
of Univci's'J. Auialytics was extremely helpjful in planing the soluta.on path and 
other jubst''ucturing items. Phase one def-Lnes the basic substructures without 
load. Phase two ocmb.lnes components and roducen thot.i to a fi"Ll pseudostructure 
but applio.'' deformations for a loading. Phase three recovers hifluenco coeffi- 
oio -.t i'latricos for the responses of individual component substructures to the 
deformations of the pscudostructure. The noxt phase is the direct transient 
analysis of the sca2.ar model of the p^cudostr*ucture, but is not caHled phase 
four. Phase four recovers the response of individual component substructures 
to the transient excitation using modified static rigid format. A post-pro- 
cessor was vjritten to scan the stress data for the times and locations of parts 
exceeding a threshold value. 
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'SCLimCN STRATEGY 


Answers need to be found tn none iTPportent quoctions ''(hen one proposes to 
uEo slfbstructur'' nc with tronoionto^ For Inctnnoo, 

V/lnt can and what enn't r ’.".■!:'’'rinted Fubntrucitiu'in^ do? 

How in the tr^rsfer of matrix data from phor.'. fwc to ir'rid ■rnv.'irt 9 
o '^coTTT.iirhcd? 

V,1n.t boo to >’t> done to the tromlont ricld format to roeccui^.o the 
mntrlco-i offnred to it bv (rubnt'‘".oturinc? 

How in dnr^vinc irt'’oduccd? 

Hov: in transient data to be roooverod ■in the componont oubstruotures? 

How arc pints to be obto^nod? 


Outline of Se’’ntion 

An vsLth any : nvostiqntlon, '’!''ow''"s to questions beget more quectlons. 

The nnswors vrill therefore not be clesr-cut until a chain is satisfied. In 
short, automated rutstruo-'-urir.c ann orgonise the oliaractoristicB of the peoudo- 
structru'Cv as it is modeled for traJ'.niontn , into its stif'^’ncso md mass matrices, 
but v/ill joLeld ncit''or a damping matrix nor load'ing. At first blush, it seemed 
feasible to reprocent this psoudontruc-^uro in sorlnr form for rigid format 9* 
Knov-flng tho rote-tned degreos of freedom in the pneud 5struc-'”.:rc, one then calls 
for an CLiuLvalont number of scalar points for transients. During transients 
tho matrices from phneo tv/o siib.structuring can bo introduced by •> DKjIP ALTER. 

All mati'ix generating modules can be by-passed and. all matrix partitioning can 
bo circuim'-ented because all SPG*s, MPG's, and OHIT's v/ore incorporated during 
the formation of the pseudostructure ati-ffnoss matrix. Taking tho response 
output from the TRD module one ran process it by influence coefficient Ltiatrd.ces 
to recover the component responses. This sounds fairly straight forv;ard, so 
one is encouraged to tangle vdth the detailed problems. 


Data Recovery 

Looking first at the data recovery problem using influence coefficent 
matrices, it helps to think in terms of n super stiffness matrix. In phase 
tvro if a unit dioplacoment were imposed in one degree of freedom vfhlXe holding 
all others to aero and this is done for each degree of freedom iii the pseudo- 
structure, it amounts to an enforced displacement in the form of a unit diago- 
nal matrix. Translating this notion to specifics means that a unit diagonal 
matrix is needed for substitution ."s the UGV matrix after module 5DR1, Now the 
SOLVE and RECOVER modules can process the results onto tho GOP (Substructure 
Operating File)® 

Appendijc A contains tho details of hov/ these ideas v/ere implemented. The 
Operations that are important are the DMI input of the unit diagonal matrix, 
the substitution of the unity m-atrix for UGV with an ECITIV, the u.sc of a SOLVF, 
command to name the pseudostructure for v/hich the solution chain is intended, 
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anfl the une of ’’ t'^' the •lolutloti latn nut on to thn "0^ in 

on orderly fashion so thnt phase three executions ca.n readily partition the data 
from the internal substruoturine book-keeping scheme. 

A rare thing was uncovered that /xlmoGt shattered this plan. A non-open- 
endedr.ess, which is quite contrary to the original design philosophy of NASTPAN, 
was encountered. Only 100 subcaaQS were provided for in any one excrution in 
lovt-’l l6. But this pneudostructure was of order 91^, To impose a unit dis- 
pDacemant in each of the 916 degroer of freedom one at a time, meant the assign- 
ment of n static subcase for enoh of the 916 enforced dof Drmations, It was s 
rare loold.ng Oase Control packet that was as I'^^rge as the Bulk Date deck. In 
discussing this unhappy event with John McDonough of Computer Sciences, he said 
that relief of this il imitation v;ns already v'orlfod out by extending the allowance 
to JOO subcases. This was small comfort in view of the need for 5 times that 
amount of relief before this job could execute. The good nev/s v/as that John 
hn.d determined, during his investigation prepcoratory to increasing the allov/anoc 
to 300, no table restriction or other Itind of overflow condition v/oui.d be con- 
fronted if a further extension were attempted. It took two tries to dilate 
both the standEird solution and the automated cubstruc tuning section to order 
1,000, This has not been generalized yet, but a scheme ic believed to be under 
consideration which wil]. uJ.low the analyst to corrmi’jjiicate his needs bo the 
03C.AR and the SHAT by either a DI/.G cr MAGTDj\M card entry giving the size of his 
non-standard subcase turay. 

The second step in data recovery is to create a set of influence coefficient 
matrices using automated substructuring phase throe. The dimension of the i,j'l 
term of the influence coefficient matrix is 


"Displacement in the i'l Unit displacement in the 

degree of freedom of per degree of freedom of the P/S," 

a Basic S/S 

where S/S means substi’ucture and P/S means pseudostructurc. Any such maitrix 
will designated where INFL represents the matrix of terms with di- 

mension u/U=l and xxx is the subscript to denote the basic S/S b;/ name, T!iis 
influence coefficient matrix will be used in a post-transient operation to per- 
form the matrix multiplication 

to obtain timp varying displacements In a component 3/S, where U(t) is a matrix 
of the p/s response displacement vectors at each of the transient outpil tires, 
and u(t) is a matrix of the response displacment vectors of component S/3, xxx, 
at corresponding transient output times. 


In computing these influence coefficient matrices during phase three, the 
number of subcases are required to match the phase two array, Tliis means that 
in the particular x-iroblera, the phase three runs for each com onor.t substructure 
contained 'JlG subcases. iico igau’. tMs plur. i.'as (il.;:o-jt shattered by a bug. 

For econony purposes the phase three runs were subnittod as restarts of the phase 
one checlEpointod r.ms, 3ac> phase one solution v;ao 3. F. 2, and each phase 
three solution vvas R. F. 1, The pliasc tljree restarts aborted. Only by the 
chanciest strokes of good fortune did Jolin McDonough happen to liave faintly 
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retnombered thut camewhcr© ho hooi’d thc.t .:oir,o diXfioultji with reetarta were over- 
come by uoint; option 9 o" t^hc flOL cai’d, i.e* SOL lt9* v^orkedl In offec'" 
option 9 avoids the conaideratior' of the optimisation featm'oa of R* P. 

This bug appeared dui'ing res tar to involving change of rigid formate only. 

Tiiiu phauo tlii’oe operation la p.ot out of tho \/oods yet, becauee a needed 
data block from tho phase one run did not get checkpointed, Axjpendix B given 
the details of hov; to i-omcdy titia defect. Appendix B also expla'no hov/ the INFL 
matri:: io equated to the UGV nvatrix imd ic then written to disc files for the 
final solution i'ceovery. 


Before talcing up the piroblenis vvitl: bi-uiisiotito, tho third and last step in 
data recovery wi]d bo completed. Aosumiiig that the matrix of p&-* .doD'tructurc 
response displacement voctoro, U(t), hue been cuccessfully written bo a disk 
file, o!;d asGuming that tho IIIHj matrices for the yevti'al basic substructures 
have boon va-ittcr. to disk files, tVic Lank of recovering basic substi'uoturo 
time varying responses in dirulaoements, streoEOS, and forces is at hand. 

These jobs will be outside the renlni of automated substructuring except that 
they will be restarted from phase one checkpointed runs. The first task to be 
pcrfox'mod is tc sot up the case control such that the vector of response dis- 
placements at each output time slice shrD.l be considered as a static aol\ition 
case. The labelling of each oabcaso the output time proves to be a great 

convenience. Novrt , the rcstai't data Ixas to be fetched in order to re-establish 
tho Internal book-ke-eping scheme oo tli„,t the OFF (Output File Froccssor) module 
can function In on ordei-ly ra,sliion, Gince Lho product QQrFLl L^(t!y_will xn'o- 
duce a matrix that con be considered the static sclutioir iiia^i’ix M, no matrix 
generatorc or matrix partitioning is needed. The first module needed niter the 
undo formed x.lot routines is GDRS so the ALTTR packet to bring in the cntrices 
and do the matrix r.aAlfclply, can sttui't just before GDR2. The outputs from SDK2 
are than delivex'cd to OFF tc satisfy the output requests in case control. The 
stress table ic output for scSiiixing by a post processor. A normal tei'iaination 
turns contrwl over to the PLOT roui-ino r.uil exits. BetailE of ho\J these tasks 
were Implamentod arc ^iven j.i'. Appendix 0. 


Transient solution 


Problems foi' transient solution begin in pliase two of oubstructui-ing. Due 
to a bug in tho commend "IITDUCE" tho macs mati'ix v-hi.ch is produced is designated 
as squai-e not symmetric. If tliis wore allowod to go un.cox'rected transient sol- 
utions v/ould be -I times more expensive tlTiari expected, because the trailers woil.d 
telegraph to tho DCOMP and FBC moduleo that unsy.:motric x'cutines vjc'.J.d have to 
be called. .1 schemo was devised, which after much otroamllning turns out to be 
disarmingly simple, using lEROK ivc change this trailer from squai-c to symmetric. 
It was no o:.;all task to discover hew simple it could be to output the poeudo- 
otruoture stiffness an'1 mass matrices. It was a mutter of discovering that 
thci'e was a module within tho substructuring lexicon that was available for ad- 
dressing explicitly but v'hich v.eis not indivrl,dual.ly featured, Tho name of the 
module is 30F0, It is incci'poratod frequently in the list;i.ngs of DI>LAT .'JjTEF 
packets for major commands, so by studying those coimnands listings it finally 
registei’ed tliat if ''they'* ceji do it then I can dc it. The two complimentary 
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:noJaloLi ‘xnd Z0FI dcacrvo to bo given individual bi3.1ing in Llic dooumonta- 

ti.a. TiL,t'.ngr if theiv vioc f;^i‘ !’'v' L'-ullor clinngo !o jyracti'y nnd fer tLo 
ou'tv'ut of ntiff'icnr. nnri men matricou arc contaiiiod ip] ondix D« 

:i profnee ti Oic t v-kr. ralntlng to 'cT'nncii.'ntL- , i t ..’o d.c'' bo '..'cTI to pon- 
do” whhi tbo nimn arot it tie v-ay lyu'.-t th-,i .,'Ulpnt ; ...nl.t for the popurlo- 
ntiMcturo 'bould he a mtrix of dirplnccniont roii; o’-ao ■•.■o' ora '..t •, no rune o n” 
f-’e .F Ir. n. rcocih''.y vo'’or:.ty and rrciQ.” t:'.:. 'utp riny bo needed. Plots 
of the pseudostruebrn’o ere coi'tainly desiroable. It Ir, tntindabory bliot the 
trojiaionto bo able to bo restarted ul a tine eai’licr than the latest tii:io of 
the preceding run, Pesttu'ting the problem v/ith an old set of initial conditions 
on u new configuration nuat include nccclei‘ation. The liiodel is to bo scaltir 
v/itli matricoE to be delivered from the jhaso two pBCudostructare. Loads ai'o to 
be applied to this scalar model. Damping nooeb- to be .inti'oducod. 

Strange to caj", the daiiipJaig problem will be discuoaoct first. T)io reason 
for ti-ia order it; to settle the question as to whether al3. nuttriac generation 
modules con be by -passed . If uniforis otructural dumping is an acceptable 
representation of the way the structure behaves, then the dtunplng mati'ix cmi be 
generated by a scalar multiplication of the stiffneca matrix, v/hich alreadj' 
exists. This was decided on. Consequently, all matrix generation modules 
could be by-paaaed and an /iLTI'D packet could be added v.hich woilcl do the muti'i;i 
multiplication and v;hatover related parameter I’^anipulLition that wou-ld be needed. 
The lulit Data of course niust contain a faiRAM cni"d for t:ho W3 frequency. I'^lea 
Hurwits of David Ttiylor NSRDC was of inestimable value i.n helping v/ith parameter 
nianipulatioiiG and other system problems. 

riov/ the opei’atiun of bringing in the stiffness Eind nuss matrices from 
plinse tvifo is simply a mttei of using INPUTTl and renaming them with an EQUIV 
statement so that ti-anrionts can proceed along the normal chain of the rigid 
formal;. The Bulk Data of course must contain an SPOINT cai-d f:ci:taining the 
number of points cqutl to the degrees of freedom Li the pscudoatructure. 


Tlie problem of load dcfi;iitlon is mrin 




, r^omo', 


ac/ilar, t';r nr.l'ure of 


nont of load luvs a separate amplifier 
xs raiUiC Lui or.tnU'1' jOkJO. IIil^ I.I-l* 

is ncccssai'y to c:msult 


Lion time hi 

VI VwcL r f'.Uoii. 


of booh“lccu^,‘.ng, Jhethor 
Ibo Ict.d As such Lh.,t each coirrpo- 
1st cry uO the TABLEDl input data 


-L X o - ' to .X i / Om ts ■ so X L 


^r.blo to determine ;/hiul‘. scalar po'nt coxu’esponds to 
a leaded component of a geometric pc '.at, Eoriunabely, the autcrmxted subs true- 
buring output items xuiticip.vte bids need very .dcoly. Th.; name of !;Iie table 
vwdeh tahulotes the correspondence botvjecn ai; u'.tcrnc,l ilogi-ee of freedom number 
iuid tLo physical point compunont is "BUl-IMARY CP rSSUDOGTRUCTUKE CCrWSCTIVITISS . ” 
Tills table is planted in res ’onso bo the a:ialyst*o selectirn of subcommand 
"OIFITUT Option 12'* dmdr.g a OCUBTIFE operation. Of course, in this case U'le 
table to bo used is blvii associated idth L'ho final CCfTIINI] operation v;/lilch pi'o~ 
duced the firal pscudo-structm’O. If however, the opcrati.cn which pi’oduced the 
final pseudostructure cunfiguration wei’e a DEDUCE operation, then one needs to 
consult 0 pcxii' of tables:, the E:'.GC set of tables and the test PSEUDOSTRUCTURE 
CUNMECTIV IVIES table. The E.^SS is printed in response to the analyst’s se3.ee- 
tion of option 5 of the subcommand OUTPUT during a IffiDUCE opei'ution. 
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Th(i problcn cf rusbax'tirG the txani:ient integration at c. t'.iiie earlier tlian 
the laat or t\ie preceding run is cornplex. The reason, for imposing this req.uire- 
mont is that damage ia expected tc occur to the stx’ucture, but the time of dam- 
age won't be knoxm until the rosultu of a previous run are examined. If a 
strecn level is found to be oxceected, the flexibility is provided to substitute 
a replacement substructure iu. a djuungod configur’ation. The reconfigured struc- 
ture will thou be rectartod at the instant at which the sbrooo level was found 
to to exceeded. This caijability bo restart at an eai'lior t'r'ie meant that a 
chenge iii the code of the TIID module had to bo made. Simaltanoously , a provi- 
sion "SI the DI4AP listings had tc be made to allow foi' the modification of two 
book-kuexjing itema. The value cf the parameter KCOL bad to be set to tell the 
TiJD rxdule at v.iiat column in the matrix of the displacomcntu, velocities, end 
arcelur.atious the data is to be futohod to represent the j.nitial conditions at 
tho time of ro^tai't. The table TOL (Transient Output List) has to be enabled 
so that additions ciui be made to the table. The chcuigcs to the code x’ncl to tho 
DMA.P were all gcnercusly provided by ?. It. T-amidi of Computer Sciences. 

The problem of initial conditions on a changed ccrifigui\ation could now be 
liandled as a rosturt at a time specified by the MCCL and a revised TGTEP cai'd 
without having to go into an exterral dcfijiitio:’ using an IC card. Such a re- 
start also px'ovides for iui'tial accelerations to be imposed as well as initial 
velocities and displacomenxta. Guch gymnastics are possible under two conditions 
First, the changed configuration must have matrices of the same oise - (order N)» 
and in the same sequence as the o-;iginal model. Gecor.dl.y, the anaH.yst has to 
be content with the approximate values of the initial conditions on the changed 
configuration being tlie same as the ti^rmlnal values of its predecessor. Under 
certain cenditione this reconfiguring w: th substx-ucturing allows an analy.st to 
get valuable information about a nonliixeur problem using linear analycis. 

Tho final transient bask arises because of tlie po.i’ticular nature of the 
UDVT iiiatrix. The CFP arranges a triplet of value:- for^every time. Consequently 
UDVT consists of U, tJ, & fl at t^» followed by ”, U, 1 U at t^ , oi‘.d so forth 
through "^ime t^y. Only tlio ni;'t"av of tine varying d5 splacements is needed for 
the ^'^ata ..'cnovcry phase; tlio'-ofore tho iicplacemc’-^ts will h ',ve to be stripped 
from UDVT by use of the DM'I’ xitulity r.'’iTvTN. Of course, a partition.-V-ng vector 
v/ill have to he defined by FMT in the HiO’- Data, to which PARTN can refer. 

Having the matrix of displnceronts only, it junt I'cmnins to use OUTFUTl to read 
them onto an exto.-nel disk 

Unfortunately, plots of the pneudostructure cannot be obtained because it 
was defined in H. F. 9 by scalars v;''ich have no goomotrj.c proportios. ITcwevnr, 
Ijlotc of each indiiriduol basic nutastructui'c can be obtained in tho data recovery 
phone. A listing rlorg v/ith ojcjilanationc of tho ALTEH packet for H. F. 9 bo 
implement these tasks arc given in Appendix 15. 


APPLICATION 


The tecliniqne of combining automated substructuring with direct transients 
w-'s appl.iod to the analysis of a ouiwent missle. Figure 1 shows a picture of 
the hardvxare. Figures 2 through 7 depict the undeformed plots of the substruc- 
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QV poou 


P^GE IS 
unM/ry 


t'ji’on, Ficure 8 ahowi. i rhart j:;lvinc the evolution of tlio r^nulyrin fro::i '-ho 
duftnitinn of tndivi'’^'il rubntructuren, the comMrin:; irto pr>eud'^r.tv>uo1nu'on, 
t’i 'ou;jh trrinciont analyoio and fi.n'illy data recovery. chart tr ’>r.un'' ly 

caUod a substructure nr.aly.tir. tree or rinihly r, tree. In tliia roro, I'TVfovor^ 
there are five rln'iea nnctead of threo. Fism’o 9 diryloyR tVie v'etatienn in 
the tren vmder different confiy;uratlonra with n tniniria''. ''.f nnnotrtinn. 


'^he tochni-'yiQ proved tc be quite succe in 

oompai’od erite favorably with, two full scale testa. 


( 1 


the nnnlyticrJ venulta 


iiiNininirr 


Natu"ally, when one ncldcveG a ce -tui.n racar.ur:.- of sucoc'-s with e task, 
the thinlnin" p’'oco.nr docs not end abruptly. It has cccuvod to 'v.o that sotne 
thinen could hove been done differently. Cue •in ia the rrea of dyuE-ndc loads. 
Allow me bo retrocrecs for- a 'nomenh. During bho atrOyGif “rent enro won token 
to control the b.-md and the don.uity of the K matrix. Banding v;aE porfoi’mcd 
nn each substructuro, Condenoationr wo"e rerfnrmed on mstriev-'r wht)”ever poo- 
olblc v.'hon they ’..ore of small order; c.^, ObUT’s '•,'oro introduced in phase 1 os 
deeply an pc 'sible '..-ithout intorferlr.^ v.’lth connoct’ion.s ra' loads; reducing ',/rn 
used in phase two at oonnection d.ntorfacer without ’'.ntorforing with Ineds. 
Cendcnr.ationo, o.ft-r leads v;ero defined in t"an.sientG^ were c moldornd but nban- 
doned, brenure the pry-off • •".anp |••>’•rde-off.': wns n.et i mined ia*'oly evident. ?ar« 


b-sc 


n 


allol condennation.n v.oald Iv'vo to be porfTmed : 
boc.-'ii.” ■' m.' .trice.-' v/ould h'ure to be delivered to T?. F. o in uncondonned form, 
then .0 ph.'se tv/o condensation vro\<ld hnvo to follow the iiv-'.trix transi’oi' to match 
the condensation that would take nlace in transients, because data rocovex’y of 
necessity i.s based upont the mtrices in ”Sn.[SPL.lGFWlNT" form from t.ransients. 
This penalty of double condensation plus the increased density may outweigh 
the rdvantage of order reduction in transnent DCObiP and FB.So A v/aj' of oraittdjig 
tlie double condensation penalty ha.s sm’facod after tlie analysi.s. The DABEA 
lopding could have been represented ac "unity" .static loading in pha.ses one 
and two v^jith condensations i-ucludlng some of the lo.aded jioints. The "static 
load matrix" could be examined after phase two .and before transients to jirrive 
at the weighting functions different thaji 1 for the DABII/. cards, Talme ampli- 
fications would have to be adjusted accordingly. 


Another post hoc idea occurred in the area of plotting. After the tran- 
sient analysis of the scalar model was conqileto, another grid point model could 
be as.seinbled with a grid point for every scalar point wherein the other 5 de- 
grees of freedom would he constrained. An ALT3F packet to admit the static 
and deformed PLOT modules and substitute the scalar IIDVT foi’ the "BDISPLACE- 
MENT" ini5ut to the VDP module would allow the plot;-, to be- obtained for the 
pseudostructui’o transient response at a small cost. 
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CONCLUDIHCf .REKu\RICS 


A tochniquG for con:binin£ automated substructure analysis witli transient 
analysis Ivrs been devised and successfully applied. The teclinique can be sum- 
mai'ized in the form of a recipe to assist an areffe^st, v/ho may want to ti’y this 
technique, as to what factors have to be taken into accourit. 


Recipe 

1. Run G/2 tlirough 02. Combine & reduce tc produce the transient P/S 
model. 

2. Run 02 SOLVE 3; RECOVER operations with /XTER packet to read ^Ij 
matrix into UGV lat^ block and supply subcases for each of the*^N 
degrees of freedom. 

3. Run 02 tc change M-ITX from sq. to sym, using MERGE* packet. Read 

TOfX S; MirfX onto external disk file usLig command and OUTPUTl 

packet. 

4. Set up transients as a scalar problam. Read Kl-fTX & I-MfX in fi'om ex- 
ternal fil,e. Supply loading data to scalar degrees of freedom. 

5. Apply /iLTEE packet to jump around ms.trix generators, to build damping 
matrix, to set ti'ansiont parameters, to partition displacements from 
UDVT, and to provide for restai't @ earlier time than last.** 

6. Run 03 as S0L 1,9* with ALTER packet to correct SDRl* and output INFL. 

7. Run 04 as S0L 1,9* vrith .\LTER packet to jump around matrix generators, 
to multiply for basic 3/S transient disijlacoment response, and output 
stress table. 


N is limited to 3^30 in ebando.rd Ic-'vc-. 17. 

* Tlrese steps will not be necess'-ry i:: ai'. operatii'ig vex-sion with these 
bugs corrected, 

** Level 17 has this feature incorpen-ated i:i bhe code. 
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Figure 3 
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APPENDIX A 


DMAP ALTER TO R.F.I. FOR 
UNITARY MATRIX LOADING WITH S/S SOLVE 

This appendix supplements the description of SUBSTRUCTURE 
Phase II Solution Strategy. The Phase II DMAP ALTER statements 
are given first, followed by explanations according to statement 
numb er . 

1. ALTER 100 

2. JUMP TGB $ 

3. ALTER 126 

4. LABEL TGB $ 

5. EQUIV PARTUGV, UGV/ ALWAYS $ 

6 . ENDALTER 

Statement Explanations: 

1. Since these ALTER statements are used in conjunction 
with S/S commands SOLVE and RECOVER, this ALTER must not inter- 
fere with the automated ALTER' s associated with these commands, 
i.e. these must be avoided 2, 4 to 5 , 9 to 22, 29 to 30, 41, 58 
to 61, 73 to 78, and 134 to 164. The ordinary preparation of 
matrices such as constraints and omits will in general be needed. 
The operations to be eschex^fed as unnecessary and costly are de- 
compositions, load assembly, FBS (Forx-jard Backward Substitution), 
and 3DRI (Stress Data P..e c every ) . Begin the ALTE?v at statement 
100 in R.F.I because constraints and omits have been completed 
by this point, and it precedes the decomposition. 


318 



2,3,4. Start the avoidance of unwanted operations with 
a JU>IP command to a label located after the recovery of UGV in 
module SDRl. 

5. The name o the unitary matrix input via DMI cards 
is PARTUGV. Since it is de.iired to have NASTRAN recognize this 
unitary matrix as the displacement solution, it must be labeled 
UGV for succeeding modules to so recognize it. There PARTUGV is 
ec, nival ent to UGV. The parameter ALWAYS was used to make this 
EQUIV "take" because it was established earlier in the automatic 
ALTER' s with a value of negative one. 
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APPENDIX B 


DMAP ALTER FOR CORRECTING SDRl AND OUTPUTTING INFLUENCE 

COEFFICIENT MATRIX 

This appendix supplements the description of SUBSTRUCTURE Phase III Solution 
Strategy. The Phase III DMAP ALTER statements are given first, followed by 
explanations according to statement number. 


1. 

ALTER 126, 126 

2. 

COND 

NOUOOV, OMIT $ 

3. 

FBS 

LOO, ,P0/U0DK/C,N,1/C,N,1/C,N,1/C,N,1 $ 

4. 

CHKPNT 

UODK $ 

5. 

EQUIV 

UODK, UOOV/ALWAYS $ 

6. 

CHKPNT 

UOOV $ 

7. 

LABEL 

NOUOOV $ 

8. 

SDRl 

USET, PG, ULV, UOOV, YS, GO, GM, PS, KFS, 



KSS,/UGV, PGG, QG/V,N,NSKIP/C,N, STATICS $ 

9. 

OUTPUT 

1, ,,,,//C,N,-l/C,N,3 $ 

10. 

OUTPUT 

1 UGV,,, ,//C,N/0/C,N,3 $ 

11. 

EXIT $ 


12. 

ENDALTER 


Statement Explanations : 

1. This is a restart of a SUBSTRUCTURE Phase 1 checkpoint. A required 
input data block UOOV for module SDRl did not get checkpointed in the 
Phase I run because SSG3 was automatically ALTERed out, so therefore is not 
available for restart and its absence would cause ari abort because UOOV is 
not allowed to be absent as an input data block if an OMIT occurred. 




Since the SDRl module in R.F. 1 is statement number 126, the 
patch for remedying this defect must be inserted cither ahead 
of or in place of 126. In this case it is inserted in place of 
126, because constraint forces QR were not checkpointed in the 
Phase I run either, necessitating the removal of QR as an input 
data block to SDRl, 

2 thru 7, These six statements are intended to provide 
for both the presence and absence of a Phase I OMIT. If there 
were a Phase I OMIT, the parameter named OMIT would be given a 
positive value when output from GP4 and the COND conditional 
jump to label NOUOOV would not take; thus causing the operation 
sequence to pass to statement 3, PBS. If there were no Phase I 
OMIT, the parameter named OMIT would have retained its nega- 
tive one default value, and the COND conditional jump to label 
NOUOOV would "take" and the PBS module would be by-passed. 

The reason that it is desired to cause this conditional 
jump is that UOOV would have existed albeit as a purged data 
block, and would satisfy the input requirements of SDRl, The 
EQUIV of data block UODK to UOOV uses parameter ALWAYS because 
ALWAYS V7as defined to be negative one by automated substruc- 
turing. Data blocks UODK and UOOV are checkpointed before and 
after the EQUIV because the 'DMAP compiler is very particular 
about tidiness in specifying the precedents and antecedents in 


EQUIV. 

3. Data block PO is required to be supplied by the 
analyst as DMi input. It is of order "A" rows and "0" colijmns 
and is null, This will create a null UOOV, The reason that UOOV 
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should be null is that no load has been put on any omitted points. 
8. The call statement for SDRl must be written without an 


input data block appearing after KSS, It is not necessary to 
remove output data block QG even though QR is absent, because 
SDRl will generate QG as purged in the absence of QR, 

Output data block UGV will nov; be the basic substructure 
influei'ice coefficient matrix. 

9,10. These modules will output the influence coefficients 
to a MSTRAN file INP3. If the JCL is written to make a disk 
file of INP3 it should be named INFLxxx for ready identifica- 
tion . 

11, An exit is taken after outputting the IRFLxxx data 
because no more processing is needed. 

Statements 4, 5, 6 could have been omitted and UODK could 
have been used in place of UOOV ar; an input data block to SDR 1, 
because USET would have sensed that UOOV was purged during execu^ 
tion of DMAP #71 of R.F, 1 in the event of no OMIT, and would 
have relieved SDRl from requiring UOOV as input to SDRl 
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APPENDIX C 


DMAP ALTER TO R.P.l FOR POST- TRANS I ENT DATA RECOVERY 

This appendix supplements the description of STATIC DATA 
RECOVERY PHASE IV. The Phase IV DMAP ALTER statements are given 
first, followed by explanations according to statement number. 


1. 

ALTER 

39, 155 

2. 

INPUTTl 

/,,,,/C,N,-l/C,N,3 $ 

3, 

INPUTTl 

/INFLxxx, , , ,/C,N,0,/G,N,3 $ 

4. 

INPUTTl 

/,.,./C,N,-l/C,N,4 $ 

5. 

INPUTTl 

/FLIT#U, , , ,/C,N,0/G,N,4 $ 

6. 

MPYAD 

INFLxxx , FLIT#U , / TRANxxx/C , N , ) / G , N , 1/ 



C,N,0/C,N,1 $ 

7. 

PARAH 

/ /C , N ,MPY/V , N , ALWAYS/C ,N , -1/C , N ,L $ 

8. 

EQUIV 

TRANxxx, UGV/ALWAYS $ 

9. 

CHIGPNT 

UGV $ 

10. 

SDR2 

CASECC, CSTM, MPT, DIT, EQEXIN, SIL, 



GPTT , EDT , BGPDT , . , UGV . EST , , / , , OUGVI 



OESl, OEFl, PUGV1/C,N, STATICS/ 



V,N,N0S0RT2=-1 $ 

11. 

SAVE 

NOSORT2 $ 

12. 

OFP 

OUGVI, ,, OEFl, OESl, //V.N, CARDNO $ 

13. 

SAVE 

CARDNO $ 

14. 

OUTPUT2 , 

,,,,//C,N,-l/C,N,ll $ 

15. 

0UTPUT2 

OESl, , , ,//C,N,0/C,N,ll $ 

16. 

0UTPUT2 , 

, , , ,//C,N,-9/C,N.ll $ 


17. ENDALTER 
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1. The operations that are needed in this step are 
fetching of restart data, forming of basic tiine varying dis- 
placement matrix, exercising of SDR2 and OFF, plotting of 
deforme’ responses, and outputting of stress data to the 
post'processor. Call statements for SDR2 and OFF in stan- 
dard form process more than is needed here, hence their input 
and output data blocks are destined to be rewritten conse- 
quently everything from before matrix generators to just before 
structure plots can be by-passed. Hence, the ALTER control 
causes a jrunp from statement 39 to statement 155, 

2,3,4, and 5. Influence coefficient matrix [CNFLxxo;] is 
read in from disc file via INP3. Transient response FLIT#U 
is read in from disc file via INP4. 

6. The matrix multiply operation [iNFLxxx] x [FLIT#u] 

=. [TRANxxx] gives the Basic Substructure time varying 
response TRANxxx. 

7,8, and 9, A control parameter ALWAYS is assigned the 
value -1 to make the succeeding EQLIV operation ''take" so that 
the Basic Substructure time varying response TRANxxx is 
equivalenced to the static solution matrix UGV. Then UGV is 
checkpointed. 

10 and 11. Basic Substructure time varying displace- 
ments, element forces and element stresses are recovered 
according to the specifications stated in Case Control. The 
output pax’ameter NOSORT 2 is saved. 

12 and 13. The displacements, forces, and stresses from 
SBR2 are formatted for printing. The output parameter CAliDNO 
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is saved. 

14, 15 and 16. The table of stresses OESl are output 
to a FORTRAN file and stored on a disc with dataset name 
STRSx 30£. This is now in a form that the post processor 
scan program can access it. 
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nPPENDIX D 


MMTX TRAILER CHANGE AND MATRIX 
TRANSFER FROM SOF TO EXTERNAL 

This appendix supplements the description of TRANSIENT SOL- 
UTION STRATEGY. The DMAP listing is given first, followed by an 
image of the partitioning vector, then the explanations according 
to statement number. 

1. BEGIN $ 

2. SOFI /K1,M1,, ,/C,N,+l/C,N,MFLT/C,N,KMTX/C,N,MMTX $ 

3. MERGE Ml, , , ,BLNKVEC,/SQRSYM/C,N,2/C,N,1/C,N,6 $ 

4. PARAM //C,N,ADD/V,N,TRALR/C,N,0/C,N,-1 $ 

5. BQUIV SQRSYM,M2/TRALR $ 

6. SOFUT //C,N,+1/C,N,MFLT/C,N, EDIT/C, N,2/V,N,ZXX/V,N,ZXK/ 

C , N , DUMPARAM/V , N , ZXX/V . N , ZXX/ V , N , ZXX/ V , N , ZXX/ 
V,N,ZXX $ 

7. S0F0 • ,M2, , , ,//C,N,+l/C,N,MFLT/C,N,MMTX $ 

8. SOFUT //C, N,+1/C,N,TOC/C,N,SOFP/C,N,0/V,N, ZXX/V, N, ZXX/ 

C , N , DUMPARAM/V , N , ZXX/V , N , ZXX/ V , N , ZXX/ V , N , ZXX/ 
V,N,ZXX $ 

9. OUTPUTl, ,,,,//C,N,-l/C,N,0 $ 

10. OUTPUTl K1,M2,,,//C,N,0/C,N,0 $ 

11. END $ 

DMI BLNKVEC 0 2 11 916 1 +PART 

DMI BLNKVEC 1 1 

2. The utility SOFI brings information form the SOF into the 
current execution. The name of the P/S is MFLT and the two data 
blocks belonging to MFLT, KMTX and MMTX are being requested and 
being given temporary temporary data set names K1 and Ml respec- 
tively. 
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3. A merge operation of the mass matrix Ml with a null matrix is 
performed by using a 916th order null partitioning vector. The 
partitioning vector is ^asy to specify as shown on the example DMI 
card, Ablank entry is given for the first row of the first col- 
umn and implied blanks for all others. The output matrix is r\amed 
SQRSYM. The 3rd parameter states that the trailer of the output 
matrix is to be symmetric. 

4,5. A parameter is defined for use in the EQUIV statement so 
that the output mass matrix can conform to the substructuring for- 
mat for naming. 

6,7,8. SOFUT is used to purge the old MFLT square mass matrix Ml 
from the SOF. SOFO is used to output the summetric mass matrix 
M2 from the current execution to the SOF record of P/S MFLT. 

SOFUT is used to get a listing of the SOF table of contents to 
verify that the mass matrix was read onto the SOF. 

9,10. The utility OUTPUTl reads K1 § M2 onto the NASTRAN file 
INPT. 
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APPENDIX E 


t)Map modifications to dirict transients 


This appendix supplements the description of DIRECT 
TRANS' fT ANALYSIS Solution Strategy, The Rigid Format 9 DMA? 
AI.TER statements are given first, followed by explanations 


according 

to statement number. 

1. 

ALTER 

2.2 

2. 

FILE 

KGGX = TAPE/KGG = TAPE/UDVT = APPEND/ 



TOL=APPEND $ 

3. 

ALTER 

30,30 

4. 

ALTER 

33,33 

5. 

PARAM 

//C,N,ADD/V,N,NOBGG=-l/C,N-l/C,N,' 

6 , 

ALTER 

34 

y. 

cond 

LBLl, NOS IMP $ 

8. 

ALTER 67 


9. 

INPUTTl 

/. , , ,/C.,N,-l/C,N,0 $ 

10. 

INPUTTl 

/Kl.Ml, , ,/C,N,0/C,N,0 $ 

11. 

EQUIV 

Kl, KGG/NOBGG/Ml, MGG/NOBGG $ 

12, 

ALTER 

102 

13. 

PURGE 

MAA/NOICGGX $ 

14. 

ALTER 

110 

15. 

JUMP 

LBL5 $ 

16. 

ALTER 

113 

17. 

PURGE 

K4AA/NOKGGX $ 

18. 

ADD 

1CAA,/K4AA/C.N, (0.03. 0.0)$ 

]9. 

CHKPNT 

1C4AA $ 

20. 

ALTER 

139 
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21. PARAIH 

22 . ALTER 

23. EQUIV 

24. EQUIV 

25. ALTER 

26. PARTN 

27. OUTPUTl 

28. OUTPUTl 

29. SDR2 

30. SDRS 

31. CHtOPNT 

32. OFP 

33 . SAVE 

34. EXIT 

35. ENDALTER 


//C,N,ADD/V,N,N0SIMP/C,N,1/C,N,0 $ 

163,163 

PPT, PDT/NOSET $ 

PDT, PD/PDEPDO $ 

167 

UDVT, USTRIP,/,,FLIT0U,/C,N,7/C,N,1/C,N,2/C,N,2/ 
C,N,2/C,N,2 $ 

,,,,//C,N,-l/C,N,4 $ 
rLIT0U,,,,//C,N,0/C,N,4 $ 

CASEXX, CSTM, MPT, DIT,,EQDYN, SILD,,, BGDT, TOL,,, 
EST,, PPT/OPPl,,,, C,N,TRANRESP $ 
0PP1,,,,,/0PP2,,,,, $ 

0PP2 $ 

0PP2, , , , ,//V,N, CARDNO $ 

CARDNO $ 

$ 


Statement Explanations: 


1,2. The TOL (TRANSIENT OUTPUT LIST) must be enabled to that it can be 


appended for any continuation of integration. Initially, when there is no 
restart, the TOL need not be appended, but for all subsequent restarts it 
needs to be appended, so it is enabled here to be ready for all runs. Instead 
of inserting a separate APPEND statement, it becomes more concise to add the 
APPEND to the existing FILE statement. Consequently, the existing file state- 
ment was removed and put back with the added APPEND statement. 
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4.5 Subsequently it is desired to incorporate a structural 
damping matrix and no other kind of damping, but to equip the 
code via a DMAP ALTER aiid not through a matrix generator module. 
Ordinarily the EMG module vjould sense the absence of viscous 
damping and automatically set the parameter NOBGG to negative 
one so that the viscous damping matrix, BGG, Xizould not be 
generated. But ixi this solution path EMG will be bypa.ssed by the 
conditional jump of statement 7 (because stiffness and mass 
matrices are input via statemeiits 9, 10, 11). Something must be 
done about parameter NOBGG so the OSCAR can do a proper job in 
providing storage space for input data blocks. NOBGG must be 
preset to -1, so the parameter entry C,N,1 must be changed to 
C,N,-1. The original PARAM statement must be replaced by this 
revised statement therefore the ALTER 33,33 was used. 

3 thru 7. All matrix data is being read in from disc and 
no matrix generator will be called upon. The compiler recognizes 
this and sets NOSIMP to -1 during execution of TAl, This is all 
very well but if the sequence of operations x^ere left unchanged 
NOSIMP x^ould cause the conditional jump of statement 30 to engage 
statement #62 next and xvould completely bypass the PARAM state- 
ment #33 in spite of all the modifications discussed above. 
Consequently, the conditional jump based on NOSIMP was caken out 
f ■ r sition 30 X'jith the ALTER 30,30, and restored after the PARAM 
executes in position 33 and needs to have a value of -tl later on. 

8 thru 11. Stiffness and mass matrices from laase II need 
to be introduced into the DMAP stream after the matl■i^: operators 
and before the matrix partitioners implying after statement 7#6,5, 
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SMA3 , and before staLemenl: #75, GSPS. Here, ALTER 67 brings it 
in jnst before GP4. 

The matrices on file FLITIil't from Phase II can be internally 
subscripted to any variety of numbers from Kl, Ml to K9, M9 and 
higher depending on the execution order of commands in Phase II. 
In anticipation of this, statements 10 and 11 use subscript 
one followed by equivalences to KGG and MGG which will always 
’’take” because parameter NOBGG was preset to -I. If K & M are 
written out with other than subscript one, the acting values 
should be v?ritten into statements 10 and 11, 

12. and 13. If there V7ere OMIT’s the SFA (Segment File 
Allocator) could not provide for MAA, because it was purged 
back in statement #28. This purge at #28 took place in this 
instance, because NOSIMP is negative as explained above. In 
order to equip SFA to provide for SMP2, YiAA has to be unpurged 
befox'e statement #103 when SMP2 goes into operation; therefore 
the unpurge is introduced at #102 with the control parameter 
NOKGGX which w^as set equal to +1 at statement #31. 

14 and 15. It is the intention to provide for the damping 
matrix K4AA by a DMAP ALTEPv. subsequently, besides which K4FF 
does not exist to act as an input data block for SKP2; therefore 
SMP2 and CHKPNT are lumped around. 

16, 17, 18 and 19. Now the damping will be generated. 

The output data block name will be K4AA, Bur this ‘.ras purged 

along with MAA at statement #28 so it also has- r:- be ur.purged; 
again the choice of control parameter is NOKGGX Srarially 
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uniform s iruciiiral clamping is proportional to stiffness so a simple 
scalar multiply will serve. The ADD statement performs 
ctA '1- CSB = c. If B=P , then C=aA. The coefficient a on A is a 
complex number, demanding that the real and imaginary parts be 
supplied. In this instance, the real part is 0.03 and the 
imaginary part is 0.0. 

20 and 21. KDEK2 is a parameter which indicates to module 
GKAD whether or not matrices ICAA and MAA are present according 
to its value of d'l or -1. Its value is computed in statement 
#140 by a logical AND operation on the values of NEGEKL aid 
NOSIMP from their most recent update according to the VPS table. 
Unless something is changed at this point NOGENL and NOSIMP 
are both " -1 so their A.ND result = -land KEDK2 will he negative. 
This negative value would indicate that KAA and MAA do not exist, 
but they do exist, because matrices were read in without benefit 
of general elements or from element matrix generators. A non- 
negative value reflects the condition of these matrices. By 
changing NOSIMP to +1, the computed value of KDEK2 is -f-l , satis- 
fying the non-negative requirement . 

23 and 24. EQUIV statement #163 has a bug in it. To 
correct this bug, the two EQUIV statements are written in place 
of the defective statement. 

26, 27, and 28, These statements provide for stripping 
velocity and acceleration vectors from the UDVT matri;;. The 
partitioning vector USTRIP is supplied by the analyst via DMI 
card data. The disx ' aceraent vectors in UDVT occur In columns 
1, 4, 7, the last output displacement. The last output 
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displnceuiinit: is cnlculniod to be; 


160 / the nvunbor of integration^ x3 / the number of vec- 


tiiiic steps selected on 
the T STKP card 


' ors per time step 
UtJ.'l/ y 


/ the out 
yselecte 


utput interval 

ted on the T STEP card 


( 3 vectors \ 
for the \ 
;jeroeth ) 
time step / 


« i'3 " 63. This value is put in field 8 of the DM1 

card indicating that the DMI vector for controlling the colunjis 
partitioning will have 63 rows The matrix of displacement 
vectors will be named FLIT #U, where y/ indicates the restart 


number starting with the zeroeth, 

29 thru 35. The call to SDR2 allowr the modules to respond 
to the Case Control Card OLOAD. Modules are limited to the input 
and output data bl 'ks necessary to process only the loads. 

Exit after OFp because no plots are requested. 
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